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A  COMPUTER  PROGRAM  TO  IMPLEMENT  THE  CHEN  METHOD  OF  DIMENSIONAL  ANALYSIS 


Introduction 

Dimensional  analysis  provides  a  method  of  reducing  the  experimentation  required  to 
analyze  a  problem  of  several  variables  by  illustrating  the  nature  of  the  relationship  between 
the  variables.  Various  methods  exist  to  determine  dimensionless  parameters,  but  the  methods 
all  have  in  common  the  following  requirements: 

1.  All  pertinent  variables  of  the  problem  must  be  included  in  the  analysis, 

2.  No  extraneous  variables  may  be  included  in  the  analysis, 

3.  Only  a  single  relationship  exists  between  the  variables  of  the  problem. 

The  procedures  for  determining  dimensional  parameters  do  not  require  knowledge  of  the 
relationships  between  the  variables  of  the  problem  but  the  results  of  the  procedures  are  often 
improved  if  some  knowledge  is  applied. 

Assuming  that  sufficient  information  about  the  problem  is  available  to  satisfy  the 
requirements  of  the  analysis,  a  relationship  between  the  variables  may  be  written  in  terms  of 
some  unknown  function  (4>)  as: 

*(qi,q=/q3, . . •  ,q«)  =  o  .  (i) 

where  the  qi's  are  the  variables  of  the  problem. 

Dimensional  analysis  provides  a  method  of  reducing  the  problem  of  Equation  1  to  the  form: 

*(tu  ,ti2  ,n3  ,  •  •  •  ,Ttm)  =  0  .  (2) 

where  the  Tti  '  s  are  dimensionless  parameters,  and  m  is  less  than  n. 


Various  methods  for  determining  dimensionless  parameters  are  available  in  the  literature, 
most  notably  the  Buckingham  Pi  method  and  the  Lord  Rayleigh  method.  These  methods  may  be 
somewhat  cumbersome  when  the  number  of  variables  becomes  large.  Chen  [1]  recently  published  a 
method  for  determining  dimensionless  parameters,  which  uses  matrix  algebra  and  is  amenable  to 
computer  solution. 


Chen  Method 

The  Chen  method  allows  for  the  development  of  dimensionless  groups  using  standard  matrix 
algebra.  The  method  is  illustrated  here  using  an  example  provided  by  Chen  [2].  Chen  considered 
the  problem  of  forced  flow  inside  a  rough  circular  pipe.  The  problem  statement  is  given  in 
terms  of  an  unknown  function  4>  of  the  important  variables  of  the  problem  by: 

$(  Ap/L,  D,  p,  d,  v,  L,  k)  =  0  .  (3) 


where: 

A  p/L  -  Pressure  drop  per  unit  length  of  pipe.  [  A  p/L]  =  MXL"2T"2 
D  -  Pipe  Diameter.  [D]  =  M°LXT° 
p  -  Fluid  dynamic  viscosity,  [p]  =  MXL*XT'X 
d  -  Fluid  density,  [d]  =  MXL'3T° 
v  -  Fluid  mean  velocity,  [v]  =  M°LXT~X 
L  -  Distance  from  entrance.  [L]  =  M°LXT° 
k  -  Pipe  roughness,  [k]  =  M°LXT° 

In  the  notation  above  expressions  of  the  form  [a]  =  MXLXT'2  should  read  "a  has  units  of 
mass  and  length  per  unit  time  squared.  Chen  showed  the  dimensional  analysis  problem  in  this 
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case  may  be  cast  in  the  following  form: 


The  terms  in  the  column  matrix  are  strictly  logarithmic  terms,  and  Equations  4  and  5  are  a 
statement  of  the  dimensions  of  the  variables  of  the  problem.  Chen  termed  the  matrix  of  Equation  4 
the  repeating  variable  matrix,  since  the  column  matrix  of  Equation  4  contains  the  repeating  terms  of 
the  Buckingham  Pi  method.  Note  that  Equation  4  has  the  form  AX  =  b,  so  that  X  =  A'^b,  where 
A-1  is  the  inverse  of  the  matrix  A.  Equation  5  has  the  form  Bx  =  C,  and  consequently 
BA-^b  =  C,  so  that: 


where  A’1  is  the  inverse  of  the  matrix  of  Equation  4  and  b  is  the  column  vector  of  Equation  4. 
Expanding  Equation  6  gives: 


Carrying  out  the  matrix  multiplication  gives: 
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Equation  8  provides  four  dimensionless  groups  as  follows: 


dv2L  |jvL 

Hi  =  TI2  = 

ApD  ApD2 


n3  =  L/D  tx4  =  k/D 


The  Chen  method  has  been  programmed  in  the  Fortran  computer  programming  language  for 
execution  under  the  Unix  operating  system.  Examples  of  the  use  of  the  computer  program  are 
provided  in  the  following  section.  A  listing  of  the  computer  program  is  provided  in  the 
Appendix. 
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Examples  of  Computer  Program  Use 

The  computer  program  was  developed  to  operate  under  the  Unix  operating  system. 
The  program  reads  from  the  standard  input  and  writes  to  the  standard  output  A  data  file 
of  the  input  data  is  usually  created,  and  input  redirection  used  to  execute  the  program. 

The  format  of  the  input  file  is  discussed  in  the  examples. 

EXAMPLE  1  --  Drag  Force  on  Sphere 

The  variables  of  interest  in  this  problem  are: 


Dimensions 

Quantity 

Symbol 

F 

M  L 

T 

Drag  Force 

D 

1 

0  C 

0 

Diameter  of  Sphere 

d 

0 

0  1 

0 

Density  of  Air 

P 

0 

1  -3 

0 

Viscosity  of  Air 

P 

1 

0  -2 

1 

Velocity  of  Air 

V 

0 

0  1 

-1 

Newton's  Proportionality 
Constant 

9c 

-1 

1  1 

-2 

The  table  above  is  a  statement  of  the  dimensions  of  the  variables  of  the  problem. 

To  carry  out  the  dimensional  analysis,  it  is  necessary  to  determine  which  of  these  variables 
will  be  the  repeating  variables.  The  selection  of  repeating  variables  is  not  completely 
arbitrary;  the  repeating  variable  matrix  must  not  be  singular  or  ill-conditioned.  The  computer 
program  tests  for  these  conditions,  however,  and  improper  selection  of  repeating  variables 
will  result  in  an  error  message  from  the  computer  program. 


5 


O  O  O  c-l  rH  CM 


Input  Data 


The  general  format  of  the  input  Tile  is: 

Number  of  variables  of  the  problem. 

Number  of  variables  in  the  repeating  variable  matrix 
Primary  quantity  number  1 
Primary  quantity  number  2 


Primary  quantity  number  n 
a  b  c  d  ...  n  Variable  name  i 

a  b  c  d  ...  n  Variable  name  j 

a  b  c  d  ...  n  Variable  name  k 


a  b  c  d  ...  n  Variable  name  m 

where  a,  b.  c,  etc.  are  the  dimensions  of  primary  quantities  in  variable  i,  j,  etc.  In  (he 
example  above,  the  primary  quantities  are  force,  mass,  length,  and  time;  the  variable  g* 
has  primary  quantities  MLT'^F1. 


For  this  problem,  the  input  file  is: 


6 

4 

force 

mass 

length 

time 


-1 . 

1 . 

1 . 

-2. 

9c 

0. 

0  . 

1 . 

-1  . 

V 

0. 

1 . 

-3  . 

0  . 

P 

0. 

0. 

1  . 

0  . 

d 

1  - 

0. 

-2  . 

1  . 

h 

1 . 

0  . 

0. 

0  . 

D 

Note  that  the  order  of  the  input  file  assumes  that  the  repeating  variables  for  this  problem  are  g*.-, 
v,  p,  and  d.  Alternate  formulations  of  the  problem  may  not  produce  the  same  results. 
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Program  Output 


number  of  variables  for  this  problem:  6 

number  of  variables  in  the  repeating  variable  matrix:  4 


REPEATING  VARIABLE  MATRIX 


force 

mass 

length 

time 

-1.0 

1.0 

1.0 

-2.0 

9-= 

.0 

.0 

1.0 

-1.0 

V 

.0 

1.0 

-3.0 

.0 

P 

.0 

.0 

1.0 

.0 

d 

INDEPENDENT  VARIABLE  MATRIX 

force 

mass 

length 

time 

1.0 

.0 

-2.0 

1.0 

P- 

1.0 

.0 

.0 

.0 

D 

INVERSE  OF  THE  REPEATING  VARIABLE  MATRIX 

-1.0000 

2.0000 

1.0000 

2.0000 

.0000 

.0000 

1.0000 

3.0000 

.0000 

.0000 

.0000 

1.0000 

.0000 

-1.0000 

.0000 

1.0000 

SOLUTION  OF  THE 

PROBLEM 

gc 

V 

P 

d 

-1.000 

1.000 

1.000 

1.000 

P 

-1.000 

2.000 

1.000 

2.000 

D 

The  drag  force  on  a  sphere  in  a  uniform  flow  field  may  expressed  as: 


vpd  Dqc 

'  v2pd2  )  =  0 

Where  the  first  dimensionless  parameter  is  the  Reynolds  Number,  and  the  second  dimensionless 
parameter  is  the  coefficient  of  drag  for  the  sphere.  The  problem  is  now  understood  to  be: 
the  coefficient  of  drag  on  a  sphere  in  a  flow  field  is  a  function  of  Reynolds  number  only. 
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EXAMPLE  2  -  Centrifugal  Pump  Performance 


Sclmittger  [3]  considered  the  problem  of  centrifugal  pump  performance;  the  variables 
of  interest  are: 


Quantity  Symbol 

Pump  Flow  Rate  Q 

Pressure  Rise  in  Pump  Ap 

Pump  Impeller  Diameter  D 

Pump  Speed  N 

Power  Input  to  Pump  E 

Density  of  Fluid  rho 

FLUID  Dynamic  Viscosity  p 


Dimensions 


M  L  T 

0  3-1 

1  -1  -2 

0  10 

0  1-1 

12-3 
1-3  0 

1  -1  -1 


Input  data 


7 

3 


mass 

length 

time 

0  . 

0  . 

-1  . 

N 

0. 

1 . 

0. 

D 

1 

3. 

0. 

rho 

1 .  - 

1  . 

-1 . 

P 

1 . 

2. 

-3. 

E 

0. 

3  . 

-1 . 

Q 

i_  _  _ 

1 . 

-2. 

A  p 
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Program  Output 


number  of  variables  for  this  problem:  7 

number  of  variables  in  the  repeating  variable  matrix:  3 


REPEATING  VARIABLE 

MATRIX 

mass 

length 

time 

.0 

.0 

-1.0 

N 

.0 

1.0 

.0 

D 

1.0 

-3.0 

.0 

rho 

INDEPENDENT 

mass 

VARIABLE  MATRIX 
length 

time 

1.0 

-1.0 

-1.0 

1.0 

2.0 

-3.0 

E 

.0 

3.0 

-1.0 

Q 

1.0 

-1.0 

-2.0 

Ap 

INVERSE  OF  THE  REPEATING  VARIABLE  MATRIX 


.0000 

3.0000 

1.0000 

.0000 

1.0000 

.0000 

-1.0000 

.0000 

.0000 

SOLUTION  OF  THE 

PROBLEM 

N 

D 

rho 

1.000 

2.000 

1.000 

M- 

3.000 

5.000 

1.000 

E 

1.000 

3.000 

.000 

Q 

2.000 

2.000 

1.000 

Ap 

Centrifugal  pumps  may  be  compared  based  on  four  dimensionless  parameters,  namely: 

Flow  Number:  Cq  =  Q/(N  D3) 

Power  Number:  CE  =  E/(rho  N3  D5) 

Pressure  Number:  Cp  =  A  p/(rho  N2  D2} 

Reynolds  Number:  R.  =  rho  D2  N/'p. 
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An  alternate  formulation  of  the  same  problem  is  presented  below.  In  this  case  the 
variables  D,  rho,  p  and  N  are  taken  to  be  the  variables  of  the  repeating  variable  matrix. 


Input  Data 


7 

3 


mass 

length 

time 

0. 

1 . 

0. 

D 

1 .  - 

3  . 

0. 

rho 

1 . 

1  . 

-1 . 

P 

0  . 

0  . 

-1 . 

N 

1 . 

2  . 

-3. 

E 

0. 

3. 

-1  . 

Q 

1 .  - 

1 . 

-2. 

10 


Program  Output 

number  of  variables  for  this  problem:  7 

number  of  variables  in  the  repeating  variable  matrix:  3 


REPEATING  VARIABLE  MATRIX 


mass 

length 

time 

.0 

1.0 

.0 

D 

1.0 

-3.0 

.0 

rho 

1.0 

-1.0 

-1.0 

M- 

INDEPENDENT 

VARIABLE  MATRIX 

mass 

length 

time 

.0 

.0 

-1.0 

N 

1.0 

2.0 

-3.0 

E 

.0 

3.0 

-1.0 

Q 

1.0 

-1.0 

-2.0 

Ap 

INVERSE  OF  THE  REPEATING  VARIABLE  MATRIX 

3.0000  1.0000  .0000 
1.0000  .0000  .0000 
2.0000  1.0000  -1.0000 


SOLUTION  OF  THE  PROBLEM 


D 

rho 

P- 

-2.000 

-1.000 

1.000 

N 

-1.000 

-2.000 

3.000 

E 

1.000 

-1.000 

1.000 

Q 

-2.000 

-1.000 

2.000 

Ap 

Here  the  dimensionless  groups  are: 

Tti  =  N  D2  rho/p 
Ha  =  E  D  rho2/p3 
Tt3  =  Q  rho/D  p 
it*  =  Ap  D2  rho/p 
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Example  3  -  Vibration  of  a  Liquid  Drop  [2] 


Taylor  [4]  presented  the  following  problem  where  he  sought  to  find  the  lowest  natural 
frequency  f)  of  a  vibrating  liquid  drop.  The  variables  of  interest  for  this  problem  are: 


Dimensions 

Quantity 

Symbol 

F 

M 

L 

T 

Natural  Frequency 

Q 

0 

0 

0 

-1 

Diameter  of  Drop 

d 

0 

0 

1 

0 

Density  of  Liquid 

P 

0 

1 

-3 

0 

Surface  Tension 

a 

i 

0 

-1 

0 

Newton's  Proportionality 
Constant 

gc 

-l 

1 

1 

-2 

Having  assembled  the  table  above,  it  is  now  necessary  only  to  execute  the  computer  program 


Input  Data 

5 

4 

force 

mass 

length 

time 


0  . 

0. 

1 . 

0  . 

d 

0. 

1 . 

-3. 

0. 

P 

-1  . 

1 . 

1 . 

-2  . 

g<= 

1 . 

0  . 

-1 . 

0. 

a 

0. 

0. 

0. 

-1  . 

fi 
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Program  Output 


number  of  variables  for  this  problem:  5 

number  of  variables  in  the  repeating  variable  matrix:  4 


REPEATING  VARIABLE 

MATRIX 

force 

mass 

length 

time 

.0 

.0 

1.0 

.0 

d 

.0 

1.0 

-3.0 

.0 

P 

-1.0 

1.0 

1.0 

-2.0 

9c 

1.0 

.0 

-1.0 

.0 

a 

INDEPENDENT  VARIABLE  MATRIX 

force  mass 

length 

time 

.0  .0 

.0 

-1.0 

INVERSE  OF  THE  REPEATING 

VARIABLE  MATRIX 

1.0000 

.0000 

.0000 

1.0000 

3.0000  1 

.0000 

.0000 

.0000 

1.0000 

.0000 

.0000 

.0000 

1.5000 

.5000 

.5000 

-.5000 

SOLUTION  OF  THE  PROBLEM 

d  p 

9c 

a 

-1.500  -.500 

.500 

.500 

i 

The  natural  frequency  varies  inversely  as  the  three-halves  power  of  the  drop  size. 


The  dimensionless  group  may  be  written  as: 


Q2d3p 

*( - -)  =  0 

g=ct 
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Summary 


The  Chen  method  of  dimensional  analysis  has  been  shown  to  be  a  substantial  improvement 
over  other  methods  of  dimensional  analysis  both  in  terms  of  the  speed  of  implementation  of 
the  method  and  in  terms  of  the  elegance  of  the  method.  A  computer  program  to  implement  the 
method  has  been  presented,  which  should  allow  for  better  design  of  experiments. 
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character*20  name(20) , dimens ion(dunmy) 

dimension  A(dummy, dummy) ,X(dummy) ,j2(dummy) , i2(dummy) ,B(20,20) ,c(20,20) 

read(*,*)M 

read(*,*)N 

write(*,3)M 

write(*,4)N 

format ( 'number  of  variables  for  this  problem:  ',13) 
format( 'number  of  variables  in  the  repeating  variable  matrix:  ',13) 
MminusN  =  M-N 
do  75  i=l  ,N 
read(*,5)dimension(i) 
format(a20) 
do  100  i =1 , N 
do  100  j=l,N 

read(*,10)(A(i , j) ,  j=l , dummy) ,name(i) 
format(dummyf5.1 ,a20) 
continue 
write(*,*) 

write(*,*) 'REPEATING  VARIABLE  MATRIX' 
write(*,6)(dimension(i),  1=1, N) 
do  105  1=1,11 

write(*,ll)(A(i , j) ,  j=l,N),name(1) 
f ormat(dunmyf 12 . 1 , 5x , a20) 
continue 
write(*,*) 
do  110  i=l,MminusN 
do  110  j=l , N 
k  =  i+N 

read(*, 10) (B(i,j),j=l, dummy) ,name(k) 
continue 
write(*,*) 

write(* ,*) ' IHDEPEKDENT  VARIABLE  MATRIX' 
write(*,6)(dimension(i),  1=1 ,N) 
format ( / , 2x , dummy ( 6x , a6 ) , / ) 
do  106  i=l .MminusN 
k  =  i+N 

write(*,ll)(B( i , j) ,  j=l,N),name(k) 
continue 

call  minvCAfdunmy^rK^, J2,I2) 
if  (k  .ne.  0)  then 

write(* ,*) 'MATRIX  ILL  CONDITIONED  OR  SINGULAR,  k  =  ',  k 
write(*,*)'E  =  ',E 
stop 
end  if 

write(*,15) 
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15  format(/, 'INVERSE  OF  THE  REPEATING  VARIABLE  MATRIX’,/) 

do  200  i =1 , dummy 

200  write(*,20)  ( A( i , j ) ,  j=l(dummy) 

20  format(dunmyfl5.4) 

do  300  i=l,MniinusN 

do  300  j=l,N 

c(i  i  j)  =  0.0 
do  250  k=l  ,N 

250  c( i , j)  =  c(  i , j)  +  B(i,k)*A(k, j) 

300  continue 

write(* ,25) 

25  format(/,‘ SOLUTION  OF  THE  PROBLEM',/) 

write(*,26)(name(i) ,  i =1 , N) 

26  format(4x, dummy (6x,a6) ,/) 
do  400  i=l,MminusN 

k  =  i+N 

400  write(*,27)  ( c ( T , j ) ,  j=l , dummy) , name(k) 

27  format(dummyfl2.3,10x,a6) 
end 

SUBROUTINE  MINV  (A,N,E,K,X, J2,I2) 

C*  MATRIX  INVERSION  ROUTINE-FORMULATED  BY  E.  G.  CLAYTON 
C 

C  A--SQUARE  ARRAY  (SINGLE  PRECISION)  CONTAINING  ORIGINAL  MATRIX 
C  N--0RDER  OF  ORIGINAL  MATRIX 
C  E--TEST  CRITERION  FOR  NEAR  ZERO  DIVISOR 
C  K--LOCATION  FOR  SINGULARITY  OR  ILL-CONDITION  INDICATOR 
C  K=0  =)  MATRIX  NONSINGULAR. 

C  K=1  =)  MATRIX  SINGULAR  (OR  ILL-CONDITIONED) 

C  X--A  WORK  VECTOR  OF  SIZE  N 
C  J2--AN  INTEGER  WORK  VECTOR  OF  SIZE  N 

C  12— AN  INTEGER  WORK  VECTOR  OF  SIZE  N 

C 

DIMENSION  A(N,N),X(N),J2(N),I2(N) 

C 

C  INITIALIZATION 

C 

K=1 

I2(l)=0 

J2(l)=0 

C 

C  BEGIN  COMPUTATION  OF  INVERSE 

C 

DO  15  L=1,N 
C 
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ooooo  ooo  oooo 


L1=L-1 

IF(L1.EQ.0)L1=1 
BIGA=-1 . 0 


LOOK  FOR  THE  ELEMENT  OF  GREATEST  ABSOLUTE  VALUE .CHOOSING 
ONE  FROM  A  ROW  AND  COLUMN  NOT  PREVIOUSLY  USED. 


DO  5  I=1,N 
DO  1  13=1 , LI 

IF ( I  .EQ.  12(13))  GO  TO  5 

1  CONTINUE 
DO  4  J=1,N 
DO  2  13=1, LI 

IF( J  .EQ.  J2(I3) )  GO  TO  4 

2  CONTINUE 
AT=ABS(A( I , J) ) 

IF( BIGA  .GE.  AT)  GO  TO  4 

BIGA=AT 

J1=J 

11=1 

4  CONTINUE 

5  CONTINUE 

TAG  THE  ROW  AND  COLUMN  FROM  WHICH  THE  ELEMENT  IS  CHOSEN. 

J2(L)=J1 
I2(L)=I1 
DIV=A(I1, Jl) 

TEST  ELEMENT  AGAINST  ZERO  CRITERION 
IF(ABS(DIV)  .LE.  E)  GO  TO  221 

PERFORM  THE  COMPUTATIONS 

OOD  =  l./DIV 
DO  7  J=1,N 
A(I1, J)=A(I1, J)*OOD 
7  CONTINUE 
A(Il,Jl)=OOD 
DO  11  I=1,N 

IF( II  .EQ.  I)  GO  TO  11 
DO  10  J=1 ,N 

IF( Jl  .EQ.  J)  GO  TO  10 
A(I,J)=A(I,J)-A(I1,J)*A(I,J1) 

10  CONTINUE 

11  CONTINUE 
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o  o  o  o 


DO  14  1=1, N 

I F ( II  .EQ.  I)  GO  TO  14 
A(I, J1)=-A(I, J1)*A(I1, Jl) 

14  CONTINUE 

15  CONTINUE 

COMPUTATION  COMPLETE  AT  THIS  POINT 
UNSCRAMBLE  THE  INVERSE 

DO  18  J=1 , N 
DO  16  I =1 , N 
11=12(1) 

J1=J2(I) 

X( J1)=A( II , J) 

16  CONTINUE 

DO  17  1=1, N 
A(I , J)=X(I) 

17  CONTINUE 

18  CONTINUE 

DO  21  1=1, N 
DO  19  J=1 , N 
Il=I2( J) 

J1=J2( J) 

X(I1)=A(I ,  ]1) 

19  CONTINUE 

DO  20  J=1 , N 
A ( I , J)=X(J) 

20  CONTINUE 

21  CONTINUE 
K=0 

221  RETURN 
END 
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